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We present a semi-automated computer-assisted method to generate and calculate diagrams in the 
disorder averaging approach to disordered 2D conductors with intrinsic spin-orbit interaction (SOI). 
As an application, we calculate the effect of the SOI on the charge conductivity for disordered 2D 
systems and rings in the presence of Rashba and Dresselhaus SOL In an infinite-size 2D system, 
anisotropic corrections to the conductivity tensor arise due to phase-coherence and the interplay of 
Rashba and Dresselhaus SOL The effect is more pronoimced in the quasi-onedimensional case, where 
the conductivity becomes anisotropic already in the presence of only one type of SOI. The anisotropy 
further increases if the time-reversal symmetry of the Hamiltonian is broken. 
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I. INTRODUCTION 

I Spin-orbit interaction (SOI) is a promising tool for ma- 
nipulating spin degrees of freedom via electric field; 
because of that, it plays an important role in various 
novel microelectronic devices.^ In 2D semiconductor mi- 
crostructures, Rashba^ and Dresselhaus^ types of SOI 
are the most important ones. In phase-coherent dif- 
fusive systems, their dominant effect on the transport 
properties is the so-called anti-localization*-^ - isotropic 
SOI-induced correction to the conductivity tensor, lead- 

' ing to the sign change of the phase-coherent correction 
to the conductivity. Rashba and Dresselhaus SOI terms 
equally and independently contribute to the weak anti- 
localization correction. The anisotropic contribution to 
the conductivity tensor is a more subtle effect (arising in 

' the next order in the weak-disorder expansion). In an 
infinite-size 2D conductor it comes from the interference 
between Rashba and Dresselhaus SOIi^ 

A pure Dresselhaus SOI affects 2D electron systems in 
the same way as pure Rashba SOI with the same am- 
plitude. However, it is known that in a system with 
mixed (Rashba-Dresselhaus) type of SOI their action is 
not independent. This becomes most pronounced in the 
special case when the amplitude of Rashba SOI is equal 

[ to the Dresselhaus one»^ In order to highlight the inter- 
ference between the Rashba and Dresselhaus SOI, it is 
convenient to consider effects, which arise entirely due 
to this interference. An example is the anisotropic con- 

' tribution to the conductivity of a diffusive (unconfined) 
2D electron gas, which is zero in case when only one 
type of SOI is present in the system.^ 

Both the (isotropic) antilocalization correction and the 
SOI-induced anisotropic correction are phase coherent 
effects; hence it is not surprising that in a fully phase co- 
herent system both of them depend singularly on SOI 
amplitudes. In the limit of small SOI the weak lo- 
calization correction diverges both in 2D and quasi- 
ID cases, while the behavior of the anisotropic correc- 
tion depends on system's geometry: in an infinite 2D- 
system with time-reversal symmetry it remains finite 



also when SOI is infinitesimal, and in a quasi-lD case 
it diverges with the vanishing SOI. Moreover, while in 
an infinite 2D disordered slab, two different types of 
SOI (Rashba and Dresselhaus) are required in order to 
make the anisotropic component of the conductivity ten- 
sor non-zero^, in a quasi-lD geometry (see Sec lXCI l the 
SOI-induced anisotropy of the conductivity tensor arises 
also for the case when the energy spectrum is isotropic 
(i.e., when only one type of SOI is present), despite the 
fact that all dimensions of a quasi-lD sample are much 
larger than the mean-free path I of an electron. Thus, 
in the phase-coherent regime the macroscopic shape 
anisotropy of the sample results in the (microscopic^^) 
anisotropy of the conductivity tensor. The SOI is still 
required, but the energy spectrum does not need to be 
anisotropic. 

According to the theorem by Vollhardt and Wolfle 
(proved in^ for the spinless case), diffuson propagator 
poles do not contribute to the conductivity if the system 
is invariant under time reversal. We extend the validity 
of this theorem to the spinful case in Appendix |Dl From 
this theorem one may expect the appearance of uncom- 
pensated diffuson divergences in systems with broken 
time reversal symmetry, which would then result in en- 
hancement of the SOI-correction to the conductivity ten- 
sor. We indeed observe such an enhancement in the 
example of a ring pierced by a magnetic flux, where the 
time-reversal symmetry is broken due to presence of the 
vector potential (see Sec. lXIIt . 

We perform our calculations using the disorder av- 
eraging diagrammatic techniquesii2iliil2ii2 In problems 
with spin, a summation over spin indices produces huge 
expressions which cannot be handled manually anymore 
in a reasonable time. We have overcome this problem 
by developing a symbolic-calculation program^* that 
(i) generates diagrams having the requested number of 
loops, (ii) calculates the Hikami-boxes, and finally (iii) 
performs the integration over the cooperon and diffu- 
son momenta. The first two stages of the program are 
imiversal, i.e., can be readily used for other calculations 
in a diagrammatic approach. The program significantly 
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facilitates the usage of the diagrammatics, especially in 
the spinful case. 

The first part of the paper is not specific to the prob- 
lem of anisotropic conductivity. In Sec. |ll| we define 
the model which we use in our calculation, then we in- 
troduce disorder-averaged Green functions (see Sec. lIII|l 
and derive the Kubo-Greenwood formula in the Keldysh 
technique (see Sec.|y|. Then we derive the loop expan- 
sion for the diagrammatic technique in Sec. I VI l and derive 
expressions for diffusons (cooperons) in Sec. IVIII The 
second part of the paper starts with Sec. IVIII] where we 
derive the incoherent SOI-correction to the conductivity 
tensor. We then proceed with the contribution from the 
weak-localization diagrams (which remains isotropic at 
zero frequency) in Sec.|lXl The results for the anisotropic 
transport in 2D and quasi-lD geometries in the presence 
of the time-reversal invariance are described in Sec. |X] for 
the case of zero frequency &) = and in Sec.|Xl]for a) + 0. 
Finally, in Sec. I XIII we give an example of how the ef- 
fect of SOI-induced anisotropy could be enhanced by 
the time-reversal symmetry breaking terms. 

For convenience we summarize some often used no- 
tations in Tab. U at the end of the paper. 



II. THE HAMILTONIAN 

Rashba and Dresselhaus spin-orbit interaction terms 
modify the Hamiltonian as follows: 



H' = 1^ + t/: + ;j'(r), K = aY^ + hY^, (1) 

where p denotes momentum operator, a and b are Rashba 
and Dresselhaus amplitudes, and !J'(r) is the disor- 
der potential created by impurities or defects randomly 
placed in the sample. We assume that Wiy) is uncorre- 
lated: 



written as 



!J'(r)U'(r') = 



InvT 



6(r-r'), 



(2) 



where t is the mean time between collisions of an elec- 
tron off impurities, v is the density of states (DoS) at 
the Fermi level, and the over-bar indicates average over 
the different disorder configurations. The Rashba SOI is 
invariant under arbitrary rotation in the (x, y)-plane: 



= z ■ [f7 X p] s (z, a, p) = aipy 



Py + ipx 



Py - ^Vx 







-02Pz 



(3) 



where a = (cti, 02, 03) is composed of Pauli matrices, and 
R^ denotes 3x3 matrix describing rotation by an angle 

cj) around the z-axis. The Dresselhaus SOI term can be 



V'jj = (z,CtT,p) = oip^ - Oipy, 

( -1 ' 
1 







(4) 



In the coordinate system (CS), rotated by an arbitrary 
angle around the z-axis with respect to the initial CS, 
the SOI part of the Hamiltonian is transformed into 

Vs(p, a) = K(R^, p, Rla) = (z, ao + bRl.CRla, p). (5) 

In case when (p - -7t/4, Vr> becomes similar to Rashba- 
SOI^^, so that the SOI term can be written as 



CR^ = 



^-10 0) 
1 



{0 0-1) 

Vs = Vsip, a) ^ {a- b)aipy - {a + b)a2Px = sA/2, 



(6) 



where the SOI-induced spectrum-splitting A and the he- 
licity (spirality-^) operators are defined as 



K=l^pl{a + bf+p]{a-bf, 
s =2 [(fl - b)pyai - (fl + b)p.vCJ2] A"\ = 1. 



(7) 



The original disorder-free Hamiltonian in the rotated CS 
may be written in the form p^/{2m) + sA/2; it possesses 
the following eigensystem: 



lp,s) 



fsAp/2 



(a+V)Px+i{a-V)Vy 
1 



Ip), s = ±l, 
(p,s|Holp,s) = + -^-r 



(8) 



(9) 



Ap = <p,s|A|p,s) = 2 ^plia + bf + p]{a - bf. (10) 

We see that the simultaneous presence of Rashba and 
Dresselhaus SOI leads to an anisotropy of the energy 
spectrum. On the other hand, ifTOll is symmetric with 
respect to the exchange a b; conversely, the same 
is true for SOI-induced corrections to the conductivity 
tensor. We note that helicity is invariant with respect to 
the time reversal: 



02slp02 = Sp S <p|s|p). 



(11) 



In the rest of the paper we perform calculations in the 
coordinate system, rotated by n /4 in the xy plane, where 
the (unperturbed) Hamiltonian is given by 



2m 2 ^ ' 



(12) 



where !J(r) = !J'(R^^^r) is the disorder potential in the ro- 
tated coordinate system. Our assumption ||2]l that (J'(r) is 
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a 6-correlated random disorder potential is inherited by 
the disorder potential in the rotated coordinate system: 

U(r)JJ(r')= 2^6(r-r'). (13) 

The Hamiltonian (12)1 defines the velocity operator v to- 
gether with the "fictitious" vector potential A: 



so that 



n m mc 

VIC 

A = — [(a + h)a2, (b - a)ai, 0] , 
e 



— + — = — m{a^ + ¥). 

2m 2 2 



(14) 



(15) 



The strength of the SOI can be characterized by dimen- 
sionless Rashba and Dresselhaus amplitudes introduced 
as follows: 

Xa = Ipfaz/h, xt = IpfbT/h, pf = ^J2m[.i, (16) 

where is the (temperature-dependent) chemical po- 
tential. Alternatively, the SOI can be characterized by 
another set of two dimensionless parameters: 

x=7^, 6 = ^, -1<6<1, (17) 

which characterizes the "total" SOI amplitude and the 
anisotropy of the energy spectrum |(T0|| correspondingly 
The choice between two parameter sets l fT6|l and l(T7)l 
becomes important when we have to expand expres- 
sions in Taylor series. On the example of the spectrum 
splitting Ap defined in flOl l we see the advantage of the 
choice l fT7|l : while the Taylor expansion of Ap in powers 
of (x, 6) is imiform, its expansion in the parameters | [l6l l 
is non-uniform: the expansion depends on the fact if one 
expands subsequently in Xa, xy or in Xb, x„. 



where ctq is the 2x2 unity matrix, and 



(20) 



From iT9[ we see that averaged GF for the Hamiltonian 
|(12)| can be obtained from the GF of the disorder-free 
Hamiltonian by substituting infinitesimal "epsilon" in 
the denominator with (2t)~^. Thus for energies close to 
ft, the averaged GFs are very different from GFs of the 
disorder-free system; the latter are strongly modified 
due to the disorder. In fact, Eq. l fT8)l is the result of the 
summation of an infinite perturbation series. 

We calculate the imiversal contribution | |37|| to the 
conductivity tensor (see Sec.|lV]below); the correspond- 
ing momentum integrals converge in the vicinity of the 
Fermi level, so that the momentum arguments of all 
Green functions are close to pf. The assumption p ^ pp 
simplifies the SOI-term Vs in the GF-expression l fT9)l . In 
the zeroth order (in powers of Hl^p/p 1) 

Vs{-p) ^ Pf [(fl + b)ai sincp - (a - b)o2 cos cp^ , (21) 

where (/j is the angular coordinate of p. This approxi- 
mation is sufficient for the calculation of the weak local- 
ization and two-loop correction in Sec. HXl and 1X1 How- 
ever, in the calculation of the zero-loop contribution (see 
Sec. lVIIl) l. higher accuracy is required: 



(22) 



: [(fl + b)ai sincp - {a - b)o2 cos cp^ 



IV. (NON)UNIVERSAL CONTRIBUTIONS TO THE 
CONDUCTIVITY 



III. AVERAGED GREEN FUNCTION IN THE 
SELF-CONSISTENT BORN APPROXIMATION 

In our calculations, we use the electron-gas model of 
the Fermi liquid^^. In the absence of SOI and applied 
electric field, disorder-averaged Green functions are ob- 
tained from the self-consistent Born approximationii: 



(18) 



where t has small (~ EH /p.) dependence on frequency 
E. The presence of the SOI changes the expression for 
Green fimction (GF) from |(18|| into 



n>0 



in 



(19) 



= o [('^0 + Sp) gf-(p) + (cTo - Sp) ^(p)] , 



We calculate the universal (i.e., independent of the de- 
tails of the energy spectrum far from the Fermi level) 
corrections to the conductivity tensor. The latter quan- 
tity is derived in linear response in the applied electric 
field (see Sec.0. In the diagrammatic approach, the SOI- 
induced correction to the conductivity can be graphically 
represented as a sum of diagrams. A contribution of an 
individual diagram is initially expressed as an integral 
in both frequency and momentum over the combination 
of Green functions and the distribution function. We 
call such an integral universal if its leading contribution 
comes from the part of the integration space, where all 
momentum and frequency arguments of GFs in the in- 
tegrand are close to the Fermi level. In the momentum 
space this means \p - Pf\ H/l; in the frequency space 
fi\E\ < T. (pf is the Fermi momentum, and T is the 
temperature in equilibrium or effective temperature in a 
non-equilibrium case.) Then the integration in momen- 
tum space can be performed assuming constant aver- 
aged DoS V and approximating J d?pl{2nfif- ^ v ^_^dE, 
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where in 2D v = m/{2nH). According to the Fermi liq- 
uid theory, only electrons with energies near the Fermi 
level behave like free electron gas so that the effect of 
the interaction between electrons can be disregarded. 
Thus only the universal corrections are expected to give 
reasonable physical results. The non-universal contri- 
butions (i) cannot be reliably calculated and (ii) cannot 
cancel universal contributions. In the diagrammatics, 
arbitrary universal contributions to Hikami boxes can 
be calculated. Unfortunately, this is not true for non- 
universal corrections: some of them can be considered 
within the diagrammatics^-, others are too complicated 
to be calculated. The impossibility to take into account 
all non-universal contributions may leave an impression 
of imperfection of the diagrammatic technique. How- 
ever, one should note that the situation is not better rn 
the non-linear cr-modelr^ where all approximations we 
use rn the diagrammatics are required as well. 



V. NON-EQUILIBRIUM KUBO FORMULA IN KELDYSH 
TECHNIQUE 

The (mean) current density in a system, characterized 
by a one-particle density matrix (DM) p is given by 



we call "normal part of the current operator" jjv: 
(p'IIdIp) =6{p-p')e^, 

(p'IInIp) -6{p-p')e^. 

When the system is perturbed by external electric field 
E = -7 the perturbation operator is given by 



(26) 



<pWlp) = 6(p-p')2;^ 















p- 


c 




mc 


P- 


c 



p--(a + a) 



6{p - p') 

,2 



(27) 



-^A{p'\Up), (p'ljolp) = -£6(p - p')A. 



Let us denote the Keldysh-contour time ordered GF as 
G. The applied electric field affects G; the first-order 
correction is given by 



5^^^G{E,E-co) = G^ 



cn 



^E—oj 



(28) 



Expressing (l28l l in a usual 2 x 2-matrix form^, we get 
perturbation expressions for GFs in ( |25t : 



j(0 = Tr[p(0|], ]=ey, 



(23) 



6^'\Gr-Ga){E,E-cv) 



eA 



where j and v denote current and velocity operators. We 
proceed with calculations in momentum representation: 

d2p r d^p' 



^[G^^),Gr"-Gi«,Gr], 



(29) 



and 



= / (H^ / ^^U^^^P^'^^^'^^'^M' 6(i)G.(£,E-.,) = -X^f|i[(fe-fe_.,)G^^,G 



'^E-co 



5H(p,r) 

(p1j|p) = 5(p-p0e 7 , H = Ho + 6V, 
dp 

where Hq is the unperturbed Hamiltonian, and the per- 
turbation term 6V describes the applied electric field, 
see iT7[ below. It is convenient to express the DM in 
terms of GFs (see §2.1 from Ref.|2l|). 

<A|p(OIA') = {f{A;t)4>{A';t)) = -flim<A'|G<(f',OIA) = 

= lim<A'| [Gk - (Gr - Ga)] {f, OIA) = - J ^e"'"' 

— (A'|[Gi,-(GR-G^)](E,E-a;)|A), (25) 

where A = (p, s). We assume that the unperturbed 
DM p^"' is stationary (though not necessary equilibrium) 
and is characterized by energy distribution function /e . 
Thus, the zero-order DM is time-independent, and the 
zero-order GFs are homogeneous in time. 

The perturbation [see | |27|| below] affects both p{t) and 
j in l|23|l . We call the correction to j "diamagnetic part 
of the current operator" jo; the imperturbed part of I 



(30) 



+ /zE_„G|i)^G^-" - hEGlv^Gl-^] , hE = l- 2/e. 
In equilibrium 



ltanh§, En>-ii 



where T is the temperature (measured in energy units). 
We split the DM into "normal" and "diamagnetic" parts 
[like we did with the current operator in |[251l1: 



■ 2 .fi 
{\\bp^{co)\\') = '-J^'-^x 



—{A'\[{hE-hE-co)Glvi,Gl-'']\A), 



(32) 



<A|6pD(w)|A') 
2 



^i:^£i<^'i[fc-.-»<5^%'5r'- ,33, 

(/ZE - l)Glv^i&'-'^'] |A). 
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Then we rewrite | |23|| in the frequency space, substituting 



j(aO = Tr [(pC) + 6pn{co) + dpoic^)) (Jn + ]d)] 



Tr I p^°^N + p^"']D + 6pN{a})]N + 1>pn{co)]d \ , 



(34) 



where we neglected non-linear (in the perturbation) 
terms. Using Eqs. l|25ll , | |26|| , and | |32|| we can calculate \ 
in the momentum representation. From 



Tr hG^^;G^] = lTr{(r,[G|]-i - {GlVh^Glvp 



= iTr([r,-,z;,]G^)=iTr 



V] 
ru- 
in 



-6,7 TrG^ 



we conclude that 

ata; = Tr [SpDjN + p^^'Sfo] 



0. 



(35) 



(36) 



Both equilibrium^^ and non-equilibriumii contribu- 
tions to the persistent current in a mesoscopic ring are 

given by Tr [p^^'jw]/ while Tr [6pN(ii')jN] is the linear re- 
sponse to the applied electric field. 

In the rest of this section we assume that A is directed 
along the jS-axis, and we measure the charge current 
in the a-direction. Then (for arbitrary energy distribu- 
tion /e) 

Oapico) = cJ^^(c^) + a^pico), a^^(w) » o^p{co), 



icoA, 



r 1 e2 r 1 (37) 

■Tr[6pN(a;)7«] = Tr [z)„G^i)^Gi-"] , 



where we assumed that the (momentum) trace is E- 
independent. (This assumption is valid for all universal 
quantities except for the Drude conductivity, see Ap- 
pendixlAll. The diamagnetic correction to the conductiv- 
ity is given by 



icvAr, 



Tr 6pD(a;);«-6pD(0);« 



In the frequency integral dE/ (In) in | [38l l, large nega- 
tive frequencies £ fi/H give important contribution to 

the result, so that the £-dependence of t in i fTSl l cannot be 
neglected; this complicates the calculation of ( |38t . The 
SOI-dependent part of l|38] | is non-universal (the mo- 
menta of GFs are not bounded in the vicinity of pf). 
Since | [38l l does not contain products of different (re- 
tarded and advanced) GFs, it cannot contain diffusons 
or cooperons; hence it is incoherent and cannot produce 
corrections to the conductivity tensor having the same 
order in SOI, as our results below [see Eqs. | [84l l, | [87t , 
and ||96]ll. In what follows, we study the universal con- 
tribution l|37|l and do not calculate I 



An attempt to use the Kubo formula l|37t for calculat- 
ing the leading (Drude) conductivity contribution leads 
to divergences. In fact, | |37t is valid only for calculating 
corrections (due to o) i= 0, SOI, interaction, etc.) to the 
main (Drude) conductivity value. See Appendix lAl for 
details of calculating the Drude conductivity. 

We derived the universal contribution to the conduc- 
tivity tensor l|37|l for a general case of non-equilibrium sta- 
tionary distribution function. The result (|37t is the same 
as the one derived for the equilibrium case^^. Thus, we 
see that corrections to the conductivity are independent 
of the distribution f imction /e . Note that this is not true 
for the leading (Drude) conductivity I IA6I I, which does 
depend on /g. 

In what follows, we always perform calculations in 
the rotated coordinate system, where the spin-orbit part 
of the Hamiltonian is given by ||6]l. In this coordinate 
system, the conductivity tensor is diagonal in all consid- 
ered geometries; its anisotropic part is proportional to 03. 
[See the discussion after (jTSjl.] We denote the isotropic 
and anisotropic parts of the conductivity tensor a with 
symbols Ois and dan: 



O = ffisffo + OanOS- 



(39) 



In an arbitrary coordinate system, the (an)isotropic prop- 
erties of a 2D symmetric tensor can be characterized by 
two non-negative scalars ffis and |aanl - isotropic and 
anisotropic amplitudes defined by 



Ois =^Tra, Iffanl = -y/Tr [(CJ - CJoCTis)^ /2] 



Tr 



(40) 



It is easy to check that both and |aanl are independent 
of the choice of the coordinate system. Finally, we give 
explicit expressions for the anisotropic part Oan of the 
conductivity in the original and rotated by 7i/4 coordi- 
nate systems: 



in the original CS o^n = Oxy, and 



in the rotated CS = 



Orx - a 



yy 



(41) 



VI. THE LOOP EXPANSION 

It is convenient to represent the different contributions 
to the averaged conductivity in graphical form as dia- 
grams. The simplest (bubble) diagram [see Fig. [2(a)l l^ 
is produced by the Kubo formula (jSTjl by substituting 
Gr and Ga with the averaged GFs Gr [given by (|T9] ll 
and Ga - [Gr]^. The bubble has neither diffuson nor 
cooperon lines. One can proceed by connecting the re- 
tarded GF Gr of a bubble with an advanced one Ga 
by a cooperon or diffuson ladder in all possible (two) 
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q-p 



FIG. 1: Two representations of the weak localization diagram, 
cf. Fig. 4.8 from Ref. llQ. We call Fig. [1(a)] "ladder representa- 
tion", and Fig. |l(b)| - "coordinate representation". 



ways. Do ing s o we obtain diagrams depicted in Fig- 
ures [T] and [2(b)] containing one cooperon and one difhi- 
son. Adding more and more difhison or cooperon lines 
in all possible ways, one obtains an infinite number of 
diagrams. In this section we describe how the most im- 
portant diagrams can be selected out of this infinity for 
further calculation. 



A. Two ways of drawing diagrams 

hi Fig. [1] the same (weak localization) diagram is 
drawn in two equivalent representations: on the Ihs, the 
cooperon is drawn in the "ladder" form, [see the Ihs of 
Fig. 3(b)| while on the rhs the "coordinate" form [wavy 
line with two ends, see the rhs of Fig. [3(b)| is used. 

The weak localization diagram is usually drawn in the 
Ihs (ladder) -form (or as topologically-equivalent "bub- 
ble with maximally anti-crossing disorder-averaging 
lines", cf. Fig. 4.8 from Ref. [T^). Below we use the 
rhs (coordinate)-form [its advantages are discussed in 
Sec. lVII Bl belowl ; Fig.[3]gives a recipe how a diagram can 
be transformed from one form into another and back. 

A diagram in the coordinate representation consists of 
Green function boxes (bubbles, triangles, squares, pen- 
tagons, etc.) connected by wavy lines (cooperons and/or 
diffusons). A vertex of a Green fimction box (GFB) may 
be occupied by a (i) observable operator, (ii) external 
field operator, or (iii) end of a cooperon and/or diffuson 
line. 



B. Loops formed by cooperons and diffusons 

There are two important momentum scales in the dis- 
order averaging technique: (i) the (characteristic) abso- 
lute value of the momentum argument in averaged GFs 
p ~ Pf, and (ii) H/I <sc pp. {I is the mean free path of an 
electron between two subsequent elastic scatterings off 
impurities.) Momentum integrals from products of GFs 
of the form 




(b) 



FIG. 2: Zero-loop diagrams: (a) the Drude bubble and (b) the 
vertex renormalization. 



usually converge within the interval pf - H/l < p < 
Pf + h/l; hence Ti/I characterizes the deviation of mo- 
mentum argument of an averaged GF JlSl l or | [T9l l from 
Pf. The assumption that "large" momentum pf is much 
larger than "small" momentum h/l is crucial for the dis- 
order averaging technique, since (pf'/S)"^ lis its main 
expansion parameter (see Sec. lVI Cl below). 

The mean scattering free path / is also a scale on which 
averaged GFs llT8|l and l(T9l decay, e.g., in 3D GR/A(r - 
r') oc exp[-|r - r'|/Z]. We can interpret this saying that 
the length of a Green function line in our diagrams is 
/. Within the disorder averaging technique we cannot 
observe effects on scales shorter than /, i.e., to say that 
the length of a GF-line is /, is almost the same, as to say 
that this length is zero; thus we can consider a GF-line 
not as a line, but as a point. 

Now if we draw some diagram in its "coordinate 
representation" (see Sec. IVI At and squeeze all Green 
fimction lines into points, the result will contain only 
cooperon or diffuson (CD) lines forming a certain num- 
ber of loops. For example, a bubble in Fig. 2(a)[ has no 
loops (since it has no CD lines which could form a loop); 
a weak localization (WL) diagram in Fig. [l(b)[ has one 
loop, and all diagrams in Figs. l4l5l have two loops. 

The number of CD-loops is equal to the number of 
independent "small" momentum variables (which are 
~ h/l), or, in other words - to the number of integrals 
over the CD-momentum variables. 



q) n (p 



C. Comparing two arbitrary diagrams 



Let us estimate two arbitrary diagrams for the same 
physical quantity. An estimate for a GFB is oc vt''~^, 
where /; is the number of Green fimction lines compos- 
ing the GFB. Every CD line has a prefactor (47tvt)~^ 
[see Sec. IVIIBI below] . We estimate the DoS as v ~ 
m/{lnhA'^~^), where A ~ h/p^, and d is the spatial di- 
mension (d = 2 or d = 3). Let us denote Li 2, Hi 2, Ci,2 the 
corresponding number of loops, GFBs, and CD-lines in 
the two considered diagrams; the quantities hij denote 
number of GF lines in the 7th GFB of the first diagram, 
and do the same for the second diagram. 

The calculation of diagrams is often much simpler 
in the diffusion approximation - i.e., assumption that 
q*l h, where q* is the characteristic momentum of a 
CD line (i.e., "small" momentum variable). [The validity 
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of the diffusion approximation in our calculation arises 
from the assumption that q*l ~ xH h.] Sometimes a 
GFB gains additional smallness of the order of q* I /H <^ 1. 
One has to calculate a GFB in order to reveal how much 
of these "extra" q*l/H it has - this is not uniquely defined 
by the number of loops. In the following estimates we 
assume q*l ~ Hin order not to mix up expansions in two 
different small parameters: q*l/h 1 and {pfl/h)~^ <sc 1. 
Then the relation between two different arbitrary dia- 
grams is estimated as 



1st diagram / U^U [njl^ 2nvT"v-^] 



Ci 



2nd diagram j ^..^ ^ J^h.^ ^nv.'..-] [^f 

(43) 

We use the fact that L, = C, - H,- + 1 for i = 1, 2, and 
hj - h2„ = 2(Ci - C2), so that 



1st diagram 
2nd diagram 



n 



{2nY \PfI 



d-l 



L1-L2 



d > 1. (44) 



As we discussed above, apart from (pfI/K)'^ » 1, there 
is an additional small expansion parameter q*l/H <^ 1; so 
the total number of expansion parameters is two. (Later, 
in Sec.|Xl]the number of small parameters is three.) The 
loop expansion predicts how large (small) an arbitrary 
diagram is only in powers of (ppZ/fi)"^. 

To conclude, the statement that "every loop brings a 



smallness (pfZ/fi)"^" is known in mesoscopics; for the 
diagrams produced by the non-linear a-model^- it is ex- 
plained in Sec. IILS.c of Ref. IH. However, we are not 
aware of earlier articles, where this statement is justi- 
fied for diagrams within the usual disorder-averaging 
diagrammatic technique; this was the reason to include 
Sec.|yi]in this text. 



VII. ZERO-FREQUENCY COOPERON AND DIFFUSON 
IN THE PRESENCE OF SOI 

Both cooperon and diffuson^'^ can be graphically rep- 
resented as a sum of "ladders"; each "ladder" is given 
by retarded Green function Gr (bold line in our draw- 
ings) connected to the advanced one Ga (bold line) with 
some number of disorder-averaging (dashed) lines. An 
elementary building block of every such ladder is made 
of one dashed line connecting Gr with Ga- Below we 
discuss a convenient way of rearranging spin indices in 
every building block of the ladder. 



A. Separating spin indices 

One can write the expression for two GF-lines con- 
nected with a disorder-averaging line in different ways; 
we will write it either as 



p[,a' <—< p',,fi' 

^ 3 

AnvT 



2nvT 



1 

^ W L [g«(Pi)^;G^(p;)],,, lGAip'2)o,GR{p2)] , 



1=0 
or as 



,'^,a'^^Z = [G.(P:)G.(P.)]„, [G.(p;)G.(p^)]^^,^, 
= t [GK(Pi)5fGl(p;)]^^^^, [Gi(p^)a,G«(p2)] , 



a=0 a=0 



(45) 



(46) 



where the identity is used: 

3 3 

26s,s,6s3s, = '^T'oT = (^«r' ' ^" ^ '^^a.v. (47) 



Note that, differently from the first lines of ||45] | and (l46l l, square brackets in their second lines contain spin and 
momenta variables belonging only to the pair of GFs from one side of the diagram (the Ihs or the rhs from the central 
disorder-averaging line). This "separation of spin indices" effectively makes the Ihs of the diagrams in | |45l l and 
in ll46l l independent from their rhs. 
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Every diagram with cooperons or diffusons contains an infinite number of elementary blocks l|45l l and/or 
Below we will always separate spin variables in them according to | |45l l or | [46l l. 



B. Defining cooperon and diifuson 

Now let us transform every elementary building block in the diffuson series according to | |45] l: 





R 






-4e 









>- 




R 


R 


R 






A 


A 






-< 


-« 





P2'/^' 



+ . . . 



/,/'=0 



where Pi - Pj = P2 - P2 arid 
1 



E ~ [1 - ^o(q , Xi^/^lq) = ^ / Tr [a„Gl(p)a,Gi-"(p - q)], (49) 

L">0 \aji ^ ^ ' 



(48) 



47ZVT 



where Trspin stands for the trace only in spin indices. So l l45l l helped us to convert the diffuson series into geometric 
series that we could sum. Analogously we use (|46t to transform the cooperon^ series: 



+ 

PZ'^' Pl'«'^ 

3 





>- 

R 






A 
>- 


;a 
— 



P2'^' P'l'"'' 





R 


R 






A 


A 
>- 





^P2.^ 



+ .. 



P2'/^' 



(50) 



= E ^P^-n [g«(Pi)^/Ga(p;)]^^^^^, [G^(p^)a,G^(p2)]^^,^, 



l,l'=0 



where pi + p'j = P2 + P2 arid 



47ZVT 



E^cl ~[1-Xc(q)]-;, ^f(q)~/(^Tr^M(P^^^^^^^^ (^D 



From |(TT) it follows that in a system with time reversal 
symmetry 



so that 



'^2Gr/^(-P)CT2 = Gr/a(p), 



X^'*(q) = X^'*(q) and Cf=D^ 



(52) 



(53) 



The series | |48t [or | [50l l1 depend on four momentum 
and four spin variables; without external four GF lines 
only Dpj_p2 [or Cp'^^pJ is left, which is a function of two 
momentum and two spin variables. 

We call the quantities Dq'' and (f^ the diffuson and the 
cooperon, respectively. 

Diagrammatically, can be drawn in two ways. 
The Ihs-diagram in Fig. |3(a)| is more similar to the se- 
ries in (|48l l, while the rhs-diagram in Fig. 3(a) stresses 





\ , ^-'P' 3 / 



(b) 



FIG. 3: Diagrams for (a) diffuson and (b) cooperon in two 
representations. On the Ihs, the diffuson (cooperon) is drawn 
as a "ladder"; on the rhs it is drawn as a wavy line. (See also 
SeclVIAjand Fig.[ll) 



the fact that the diffuson without external four GF lines 
has only two ends. 



This rhs-diagram in Fig. 3(a) reflects better the spatial 
structure of the diffuson in the coordinate space. From 
Eq. | |13] l one can see that the distance between points 1 
and 2 in Fig.|3]is zero, and the same is true for the points 
3 and 4. We merged these identical points in the rhs- 
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diagram in Fig. 3(a) so that 1=2 and 3=4. A diffuson 
at frequency co in coordinate space decays on the scale 
min(|L„|,Ls) » /: 



L^, = II yflkJ^, Vi = (1 + 0/ V2, Ls = II VZ^, 
aJT <^H, X => min(|L„|, Ls) » /. 



(54) 



Thus, we see that the distance (in coordinate space) be- 
tween points 1 and 2 (or 3 and 4) in Fig. 3(a)| is (within 
our accuracy) infinitesimal, and this fact is graphically 
reflected by merging points 1 and 2 together (as well as 



points 3 and 4) in the rhs of Fig. 3(a) 



A similar reasoning is valid for the cooperon as well, 
see Fig. |3(b)[ 



C. Explicit 2D-expressions for q = 

The diffuson at zero momentum, q = 0, can be calcu- 
lated without assuming that the SOI amplitudes Xa,b are 
small. Using the fact that 

GKp) + Gj(-p) = Gk(p) + Gr(-p), 
one obtains a sum rule: 



where are components of the zero-momentum dif- 
fuson given by ( I57I56I I. 

Since the vertex renormalization of this type occurs 
in every diagram for the conductivity, we take it every- 
where into account by substituting the velocity operator 
v„ with its renormalized value Va = Palm', the only excep- 
tion is the zero-loop contribution (calculated in Sec. lVIIl |l, 
which is represented by the diagram in Fig. |2(a)| with only 
one velocity vertex being renormalized. 



D. Explicit 2D-expressions for q 

Because of the SOI, most components of the diffuson 
do not have a pole at zero momentum and frequency 
even if the dephasing effects are neglected; the diffuson 
gains a non-zero "mass", see | [57t . In case of pure-Rashba 
or pure-Dresselhaus SOI, this "mass" is quadratic in the 
SOI amplitude. This still remains true in case when one 
(Rashba or Dresselhaus) SOI amplitude is much smaller 
than the other one, so that the SOI-induced anisotropy 
of the energy spectrum iflTl 6 is a small parameter. Con- 
sequently, the integrals in diffuson momenta k and q 
converge on the scale of q < xhll, and it is convenient to 
introduce dimensionless variables 



■00 



^33 



(55) 



For q - 0, Xq is a diagonal 4x4 matrix with elements 



voo 



1 + (x„ + xi,y +K' " K' 
X= yJ[l + {x„+Xbr][l + {Xa-Xtn 



The components of the diffuson at zero momentum are 
given by the diagonal matrix 



^ IniT ^ 
47tvtDo = -^^Do 



= diag 



4a + 



l+K 



l + K 



K 



{Xa + Xbf {Xa-Xbf'K-l 



(57) 



where the electron (orbital) dephasing length (due 
to inelastic scattering) serves as a cut-off for the infinity, 
and the first element l?^lP' does not contribute to physical 

quantities when the Vollhardt-Wolfle theorem holds (see 
Appendix IDIi. The sum of the two diagrams in Fig. |2]is 
equal to the diagram in Fig. 2(a) with one velocity vertex 
renormalized: 



>> 




(58) 



which corresponds to the expression 

3 



m ' 



(59) 



K = IIlIxH, Q = Iqlxh, 



(60) 



sothatX < landQ < 1. Like in Sec. lVII D[ the calculation 
shows that only the (1, l)-minor (i.e., 3x3 matrix block) 
of the 4x4 diffuson matrix is affected by the SOI. In other 
words, the upper line and the left column of the diffuson 
matrix are independent of SOI: 



D"" = ^ i , 

lm% fq^llh - icoT 



(61) 



The element D^" gives no contribution to the conduc- 
tivity in systems with time-reversal symmetry (see Ap- 
pendix |D|; it becomes important, when this symmetry 
is broken, see Sec. IXIII In the rest of this Section we 
reduce 4x4 matrices of Xq and Dq to corresponding 
(l,l)-minors. 

To simplify the calculation, we further assume that 
x 1 (diffusio n ap proximation). Then the diffuson is 
obtained using (|49l l with Xu given by 



Xr 



1 



^2|y(0) _ 



where 5 is defined in i flTl l, and 



(0) 



2 2 



(1,0) 







1 -2iQ,] 
1 -2zQ, 
[2iQ, liQy 2 

iQ. ^ 

1 -iQ, 



-iQx iQy 



(62) 



(63) 



(64) 



10 



Note that the above expression for Xd is Hermitian and 
obeys the sum rule (55[ . 

The resulting diffuson matrix Dq has a denominator 
oc (det Yq')", where n > is an integer, and 



SdetyjJ) =2 + Q2 + q6 



where we used | |2T] |. The rest of l(66|l equals 



\ mc J m 
2 i—i \mcl 



n 

2[.IT 



(68) 



The expression | [65l l is independent of the direction of 
Q, so the same is true for the denominators of all dif- 
fuson components. Consequently, an arbitrary diagram 
containing a diffuson line with a non-zero momentum 
(e.g., rhs of Fig. ID has denominators [consisting of 



powers of detY^', detyj^;^^, and detYg'], which are 
invariant with respect to two "mirror reflections": (i) 

(K, -X,,Q. ^ -Qx), (ii) (^Cy ^ -Ky,Qy ^ -Qy), 

and (iii) (Kt <^ ^y,Qx ^ Qy)- [Note that the original 
Hamiltonian ((l), l(l2) l does not possess any of these sym- 
metries.] These symmetries are used in the programM 
for reducing the size of the integrands. 

Finally, we note that the easiest way to obtain the re- 
sults of this section is to utilize the computer programii. 



AO) 



/•(O)l 



VIII. THE ZERO-LOOP CONTRIBUTION 

In the zero-loop approximation (ZLA), the calculation 
can be performed without assuming that SOI amplitudes 
are small, <SC 1 (i.e., without assuming the validity of 
the diffusion approximation). Only two diagrams (see 
Fig. IJll having zero loops contribute to the ZLA. As we 
discussed in Sec. I VII CI their sum is equal to the diagram 
in Fig. 2(a) with one velocity vertex substituted by its 
renormalized value i 



(0) X ^ T 

»2 



~ h J (271, 



h 



I'rttjR — CjA - CTo 



m 



2|f^,f] (66) 



mc I m 



where g^^^ are defined in (|20l , and we subtracted the SOI- 
independent Drude conductivity od jAS} . The second 
line of | [66l l gives 



d^P PaPfl 



= 2^ [{xl + xl)ao + XaXftas] , 



(67) 



We see that the anisotropic terms in (|67|l and l|68] l cancel 
each other so that the charge conductivity l l66l l is propor- 
tional to the unity tensor—. Thus, within the ZLA Pan = 
and 



IS ^ 



h ijux/fi h 



X o[{jjT/ny 



(69) 



which is confirmed by the computer alge- 
bra calculation^- for the limiting case when 
2XflX;, + 1. 

The absence of CD-loops in the ZLA-diagrams in Fig.|2] 
means that the ZLA-contribution neglects interference 
between electrons. Thus, the result | |69|| is valid also for 
the phase-incoherent system (e.g., at high temperatures). 
The ZLA is local - that is, independent of the macro- 
scopic (on scales » /) geometrical details of the sample, 
being the same in 2D and quasi-lD cases. Since the ZLA- 
diagrams contain no crossings of disorder-averaging 
lines, their contribution coincides with the results of the 
Boltzmann equation approach, see the discussion in §9.6 
fromRef.m 

Note that at finite frequency oj there are non-zero 
anisotropic corrections to the conductivity tensor. We 
do not present them in the main text; see^^ for details. 

According to the loop expansion (see Sec. IVIjl , dia- 
grams having one (weak localization) and two loops 
may produce the contribution to the conductivity of the 
same order, or even larger, than l|69)l . We calculate con- 
tributions coming from these diagrams in the following 
sections. 



IX. THE WEAK LOCALIZATION CONTRIBUTION 

The weak localization contribution is provided by the 
diagram in Fig. |l(b)| with the renormalized (dashed) 
Hikami box 




(70) 



where both vertices are renormalized according to 
and (|59] |. The contribution of the diagrams in l ITDl l can be 
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written in the form 



where we used l l53l l and 



i.i=Q 

■ C , 
41 y-^ 



(71) 



2-1 .max 



in)^ spin 



CT,GR(k - p2)^^^-^GA(k - p2)a|G^(p2) 



(72) 



2 - 2 . max 



B«/^(k) = r -1- 

^ ^ J (271- 



«)^ spin 



^G«(p)a^G^(k-p)- 



:G^(k-p)a[GA(p) 



m 



(73) 



2-3 .max 



(74) 



The expressions ((72] |, l(73] |, l(74l l are generated by our progran^i^ and taken from the (automatically created) files 
2-1 .max, 2-2 .max, and 2-3 .max. 

I 



In the absence of orbital dephasing and at zero fre- 
quency, the isotropic part of ((TOl l diverges reproduc- 
ing the well-known result for the weak antilocalization 
correction^i^i'S^. 

The anisotropic part of ([70t, l(7T1l converges. Its leading 
(in the SOI) contribution is given by 



h 



(75) 



where we assumed that IXgXt <«c + 1, and 



cu _ 

'-'20 ~ 



(13171 + 262arcctg V7)/ V7- 88 - 7 log 2 



2247T 



0.14 
(76) 

These results are obtained ir^, where the expression 
for the renormalized Hikami box (together with other 
details of calculation) can be found. 

The results ( |69l l, | [75] l manifest the general rule | |39l l: 
the disorder-averaged conductivity tensor is diagonal in 
the considered (rotated by 7c/4) basis, where the SOI 
is given by l|6ll. To demonstrate this rule, consider an 
arbitrary diagram produced by averaging the Kubo for- 
mula | |37t for the off-diagonal conductivity element Oxy- 
Let us change the sign of px and of 02 everywhere in the 
expression for the diagram. The identity l|47l l remains 
valid if the sign of any Pauli matrix is changed, so that 
expressions for diffusons and cooperons will not change, 
as well as Hikami boxes, except for the Hikami box with 



the vertex Vx, which will change sign. Thus the total 
expression will change its sign; on the other hand, since 
our transformation is only the change of variables (over 
which the Tr is taken), the expression must remain in- 
variant. So we conclude that every diagram for Oxy is 
zero, and the disorder-averaged conductivity tensor is 
diagonal. 



X. TWO-LOOP CONTRIBUTION 
A. Expansion in SOI 

In the remainder of the paper we assume that the SOI 
parameters l fT7|l are small, x 1, |6| <s; 1 (we also as- 
sumed this in Sec. HXl l, and expand contributions to the 
conductivity in powers of x and 6. According to Ref. 0, 
in an infinite 2DEG, 



ayy(-fl, b) =Oxx{a, b) = axx{-a, -b) 
or ayJx,-d) = axx{x,6). 



(77) 



Then the anisotropic part of the conductivity tensor is 
given by the following expansion: 



CTan = 



m,n,r>0 
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FIG. 4: Three relevant two-loop diagrams, which contribute to 
SJq. See RefBfor more details. Each diagram contains small- 
momentum singularities, which mutually cancel each other in 
accordance with the theorem from AppendixiDl 






FIG. 5: Six irrelevant two loop diagrams, which do not con- 



tribute to Sjg. 



Physical quantities should depend only on even powers 
of SOI amplitudes; thus, we expect that V«, r S|„„ = for 
arbitrary odd m; our calculations confirm this statement 
for several values of n and r. The two loop diagrams 
can only affect terms with r > 1 in l(78] l: the calculation 
shows that SJq ^ 0. Below we calculate cr an (i) for a 2D 
case (an infinite film - see Sec. lXBb , (ii) in quasi-lD case 
(infinite wire - see Sec. IXCL and (iii) for the quasi-lD 
ring pierced by magnetic flux (see Sec. lXIII l. 

From Eqs. ([ZZj, GS), and gill we conclude that the 
anisotropic contribution aan to the conductivity tensor 
[see the definition in l|40l l1 can be extracted from the odd 
in 6 part of Oxx- 



Oxxix, 6) - Oxxix, -6) 



(79) 



B. The 2D case 



In total, there are nine two-loop diagrams: three are 
shown in Fig. U) and the other six in Fig. |5l The cal- 
culation in^* shows that the leading contribution Sjp is 
produced by three diagrams in Fig. ID Each relevant di- 
agram consists of two Hikami boxes (dashed) and three 
diffuson and/or cooperon lines. For Ref .[^ we calculated 




FIG. 6: The sum of these nine diagrams with non-dashed 
Hikami boxes is equal to the diagram in Fig. |4(c)| 



the Hikami boxes manually (that is, we programmed 
expressions for them ourselves and then computer evalu- 
ated them). In order to facilitate the calculation, we made 
a variable change, which allowed us to express the sum 
of three diagrams in Fig. H] as the diagram in Fig. |4(a)| 
with the renormalized upper Hikami box (HB). We do 
not use this trick in this paper because (i) it works only 
for systems with time-reversal invariance and (ii) such 
tricks became useless after we modified the program^-^, 
which now generates programs for calculating Hikami 
boxes of arbitrary diagrams. 

Let us now describe how Hikami boxes (HBs) for di- 
agrams in Fig. m are calculated. Every diagram in Fig.|4] 
contains two dashed HBs, and every dashed HB is given 
by the sum of three non-dashed HBs, like in Eq. jTOl l. 
Thus, every diagram in Fig.|4]is a sum of nine diagrams 
with non-dashed HBs. For example, the diagram |4(c)| 
can be expanded into a sum of nine diagrams in Fig. |6| 
Let us take the diagram|6]i) as an example. It is equal to 




(80) 



where we used < |53t . The quantity H^J'^^^^ is the product of 
two HBs; it is calculated in the (automatically generated) 
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filei^ batch . me/ 5 - 1 . max: 

3 5, . 



i,i=0 



p spin ' 



'/3/1S2 



Tr Tr 

p spin 



KGli-p + k2 - ki)aJ^GA(p + ki)] , 
-A(-p)(-|) 
xGR(-pXG^(p + ki)], 



(81) 



Tr Tr 

P spin L m 



G^(p)a}^G^(p-k2 + ki) 



xa/,GR(-p + k2)a;G^(p)]. 

These expressions could be simplified, but we on pur- 
pose wrote them in the same form, as they have been 
generated in the file batch, me/5-1, max in order to make 
the comparison easier for an interested reader. 

The sum of the three diagrams with dashed HBs in 
Fig-Hlis equal to the sum of 27 diagrams with non-dashed 
HBs. These HBs are calculated in 27 automatically gen- 
erated files batch.me/5-1 .max. . . batch. me/5 -27 .max 

After the HB-calculation, the integration over the dif- 
fuson momenta [k and q in ||80] |1 is performed. While 
the calculation of Hikami boxes is done fully automati- 
cally, integration over diffuson/cooperon momenta must 
be manually programmed, differently for different prob- 
lems. The technical details are discussed in Appendix ICl 
The resulting anisotropic contribution to the conductiv- 
ity is given by | [78l l for r - 2 and m - n = 0: 



(82) 



where 5 is defined in ((T7|l and the coefficient Sjg is 
given by 



in 2D 



Sjo = -5.6 X 10"^ 



(83) 



The expression ( |82t has a non-obvious property con- 
nected with the symmetry of the energy spectrum dlOt 
with respect to the substitution a b: the two limiting 
cases Xfl x;, and x„ » x;, are described by the same 
expression Using Eqs. ^ and ([ZZl, we rewrite ( 
in the form 



00 



xl + In pfl' 



2xaXi, <^ xl + xl ^ 1. (84) 



The contribution | [84|l is non-analytic for small x: an 
infinitesimal SOI coupling results in a finite correction 
a*^^ to the conductivity tensor. We note that a similar 
non-analyticity occurs also in the weak-antilocalization 
problem: if one neglects the dephasing (assuming infi- 
nite dephasing length L^j), an infinitesimal SOI reverts 
the sign of the weak-localization correction, switching 



it to the antilocalization regime. Similarly to the weak- 
antilocalization problem, the non-analyticity in ( |83t can 
be smeared by introducing the finite dephasing rate 
> 0. A convenient way to do this is to consider 

the response for the finite-frequency electric field; see 
Secim 



C. The quasi-lD case 

Consider a long wire whose cross section L± is much 
larger than the mean scattering length /, but with the cor- 
responding Thouless energy Ec± = filvj:/2L^ being large 
compared with the SOI-induced spectrum splitting (|7)|: 



Ec±T/fi » = max 



(85) 



When | |85l l holds, Hikami boxes remain two- 
dimensional, while diffusons and cooperons become 
one-dimensional (cf . §7.4.1 from Ref. lTo|) . In other words, 
in dHB still Trp = j d^p / (27i?j)^ , but in 1I8OI1 the momentum 
integration becomes one-dimensional. 



J (27tS)2 J (27lS)2 ^ Ll j_ (27lfi)2 J_, 



dqx 



and one should set ku = 



{InKf-' 
(86) 

= in the integrand. 



From (|86] | one can estimate that the quasi-lD SOI- 
induced contribution to the conductivity is /^/(xLj_)^ » 1 
times larger than the 2D one. 

A quasi-lD sample is macroscopically anisotropic 
(i.e., does not possess rotational symmetry), but this 
anisotropy becomes relevant only on the scales much 
larger than I. Thus, the Eq. ((77t , as well as the claim 
in Ref . d that "the conductivity tensor is isotropic when 
6 = 0" are valid in the ZLA and for the WL-diagram 
(since the distance between all vertices in corresponding 
diagrams is of the order of Z), but may not be valid for the 
contribution of the two-loop diagrams in Fig. ID In fact, 
we obtain the following anisotropic contribution for the 
conductivity tensor in the quasi-lD case (in the rotated 
CS): 



h pfI x^L\ W 6.7, 



-0.39 



+ 6 



-852 
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which leads to 



.(2) - 



n P_ 

h Pfl xP-L\ 



(3.5 + 4336) . 



(87) 



(88) 



Thus, in the quasi-lD case the singularity in the SOI- 
correction to the conductivity tensor is more pro- 
nounced: (i) it occurs even when the (averaged) energy 
spectrum is isotropic and (ii) it diverges at vanishing SOI 
when orbital dephasing effects are neglected. 
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XI. THE FINITE FREQUENCY CASE 

In this Section we will see that the anisotropic part of 
the conductivity tensor becomes an analytic function of 
the SOI amplitudes {Xa,Xj,) at finite frequency co + 0; we 
assume that the frequency is large enough compared to 
the strength of the SOI so that |a)T| » x^. 

Here we have to expand corrections to the conduc- 
tivity tensor not only in SOI, but in powers of \cox\ as 
well. Similarly to the zero-frequency case we choose 
three expansion parameters, Xc, xj, and W, defined in 
the following way: 



? « 1, 



(89) 



Xc = V-2f6JT, Xi 

W = lXcXlx\, \Xc\, \Xi\, |W| 1. 

Differently from the zero-frequency case (see Sec. lVIIll , 
the effect of the SOI on the diffuson can be calculated 
perturbatively in case when \cox\ » x^. The leading 
contribution is equal to the diffuson in the absence of 
SOI: 



^ n 1 



6,1/3, a,^ = 0...3. 



(90) 



The derivation of SOI-corrections to (|90] | is straightfor- 
ward, but lengthy, so we do not present it in the text of 
the article; see the program^* for more details. 

Comparing ll90l l with expressions for the diffuson at 
a; = (see Sees. IVIICI and IVIID| l we see that the dif- 
fusons for \cot\ x^ and \cl>t\ » x^ are very different. 
Consequently, while at 6j = the vertex renormalization 
cancels the anomalous part of the velocity operator, this 
is no longer the case for co i= and at large frequencies 
the effect of the vertex renormalization is negligible. 

Despite that, the calculation of the finite frequency case 
is similar to the one for co - 0; differences come from the 
fact that now we have three expansion parameters 
instead of two ((T7|l for a) = 0. We obtainii 



0^' = 

'-'an 



o -2ia)T-2x„Xi 1 
-2 ■ 0.25 ■ — r -, 

{x; + xl- HojtY 2n pfI 

2x„Xi, <^ X^ + X^ <^ COT ^ 1. 



(91) 



This finite-frequency result obtained first ir^^canbe inter- 
preted in terms of dephasing; substituting -icvz T/zcp, 
we obtain 



.(2) _ 



5.6x10" 



T- + T+ 2n jUt' 
■(---) 2^^' 



T0 T4 



(92) 



where the Dyakonov-Perel' relaxation times are defined 
as2^2T/T± = (Xfl + Xbf. 
XII. THE QUASI-ID RING PIERCED BY MAGNETIC 
FLUX 

The simplest way of breaking the time reversal invari- 
ance of the Hamiltonian is considering a constant vector 



potential field A. It arises in a quasi-onedimensional 
ring pierced by a magnetic flux. (The magnetic flux can 
be described by subtracting eA/c from momentum argu- 
ments of GFs.) In a ring geometry, the SOI-interaction 
of the type lO cannot be provided by usual Rashba and 
Dresselhaus mechanisms. However, such SOI is not for- 
bidden and thus can occur due to different reasons, e.g., 
like in InAs nanowires^*. 

Let us assume that A is directed along the ring's 
circumference so that A = = A^. The presence 
of non-zero vector potential breaks the identities l|52)l 
and l|53] | together with the generalized Vollhardt-Wolfle 
theorem from Appendix O thus making the anisotropy 
effect more pronounced because of the small-momentum 
diffuson-singularities, which now remain uncompen- 
sated. 

Eq. | |53|| now changes into 



c; 



D' 



q-2eA/i:' 



(93) 



The calculation of the HBs is the same as in Sec. IX Bl and 
Sec. IX CI Like in Sec. lXCI we have to perform summation 
over two diffuson variables, see | [86l l. For the diagrams 
which contain no cooperons, this summation is performed 
in the same way, like for the infinite quasi-lD wire, see 
Eq. | [86l l. The diagrams with cooperons become different: 
in every such diagram, eAy/c is always subtracted from 
one (out of two) cooperon momentum. The summation 
rule becomes then different from I 



J {2nrirJ {2nrir j_ (27i?z)2 27iL|| ^ ' 



dfcll 



with L|| being the circumference of the ring, and Lj_ its 
cross section. The summation is performed over all in- 
teger n and it cannot be approximated with integration; 
one can use the Poisson summation formula instead: 



(95) 



where the leading (A-dependent) contribution comes 
from the terms with n = ±1, and we will neglect con- 
tribution of terms with \n\ > 1. Keeping only divergent 
terms (i.e., terms, containing massless cooperon/diffuson 
matrix elements), and assuming that xL^ » I we obtairti^ 
the flux-dependent correction to the conductivity tensor: 



COs( 27I-— 



(96) 



with the coefficients 
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E; = fl,e ^'1 + exp 



-y T^V^^ {^'COS yV^^+l +C,sin ^ ^/^^^TT | , 



2 

2L2 + 274L|| - 219 



-/ «i = -2Lm - Lii - 1, 



128 



7124 V7L|| + 4513 a/2V2- 1 V7 - 2965 Jl^f^+l ^^"^^ 
bo = = = -3L|| - 2.1, fci = -1.9 ■ W^Lf. - 23U - 6.7L + 4.2, 



Co 



1792 V7 

28L|| - 4513 A/2V2 + I V7 - 2965 
1792 V7 



6 ■ 10"%! +4, ci = -8.6 ■ lO-^Lfl + 4.5L2 - 0.46L|| - 0.87, 
I 



where Ly = xLy/Z. Like in Sec. |Xl]we assumed that the 
divergence at small SOI (x — > 0) is regularized by the or- 
bital dephasing. In Eqs. l|97|l we wrote numerical values 
of bi and Ci instead of their analytic expressions in order 
to save space. The analytic expressions for bi and C\ can 
be found 

The result (|96l l has the same order of magnitude in 
the loop expansion [in powers of {pfll%)~^\, as (i) the 
infinite-plane result ( |82t and | |83t , as well as the (ii) quasi- 
ID result ||88] |. However, out of all three considered 
geometries it is the most sensitive one with respect to 
orbital dephasing. In fact, in the coherent limit, Lff, 00, 
it diverges as oc Lp, for finite SOI-amplitude x, and as 
oc in the limit x « //L^ — > 0. 



XIII. CONCLUSIONS 

We presented symbolic program^^ for generating, sort- 
ing, and calculating diagrams in the disorder-averaging 
diagrammatic technique. This program strongly facil- 
itates analytical calculations, allowing one to calculate 
subtle effects due to spin-orbit interaction which were 
virtually inaccessible before due to the large number of 
integrals to be calculated. The possibility to automa- 
tize the calculation improves the usefulness of the dia- 
grammatic approach, especially also in comparison to 
the non-linear 0-model^- , as a tool for studying disor- 
dered systems. 

Using this program, we studied anisotropic correc- 



tions to the conductivity tensor due to the spin-orbit 
interaction (SOI). The arising anisotropy is a phase- 
coherence effect; therefore it strongly depends on the 
geometry of the sample. In the quasi-lD wire the 
anisotropic correction is larger than in an infinite 2D- 
plane. Moreover, while in 2D-case the effect arises due to 
the anisotropy of the energy spectrum induced by the in- 
terference between Rashba and Dresselhaus types of SOI, 
in the quasi-lD case the conductivity is anisotropic even 
in the presence of only one type of SOI (Rashba or Dres- 
selhaus), that is, when the energy spectrum is isotropic. 
The (microscopic) anisotropy of the conductivity tensor 
arises due to the macroscopic (shape) anisotropy of the 
sample (on the scale much larger than the mean free 
path /). 

We also studied the case when the time-reversal sym- 
metry of the system is broken (a ring pierced by a mag- 
netic flux); then the anisotropy of the conductivity tensor 
becomes more sensitive to orbital dephasing effects due 
to the uncompensated small-momentum divergences in 
the integration over the diffuson momentum. 

In all the considered geometries, the effect is non- 
analytical in the amplitude of the spin-orbit interaction, 
if the orbital dephasing effects are not taken into account. 
Once the dephasing effects are considered, the conduc- 
tivity becomes an analytical fimction of the spin-orbit 
amplitudes. 

We are grateful to M. Duckheim and D. Maslov for 
helpful discussions. We acknowledge financial support 
from the Swiss NSF and the NCCR Nanoscience. 



APPENDIX A: INHOMOGENEOUS DISTRIBUTION OF ELECTRIC FIELD IN A HOMOGENEOUS WIRE 

The Kubo formula (|37l l is not valid for the calculation of the leading (Drude) contribution to the conductivity. From 
Sec.|V]one obtains: 

— f r-r ie^AjM C"' dE r r ri 

I^Mn] = J_ ^/^Tr [Plgfg^l a, 0, (Al) 

where the main contribution to Tr[^pagj?^aj is given by Tr [^p^^gjrg^j. (Note that we do not take into account spin 
degree of freedom in this section and use expressions i fTHI l for Green functions; consequently, Tr operators here do 
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not contain trace over spin.) 

Only values of E close to the Fermi level (E - 0) contribute to the integral in IIAlll , so that 



= IbfiT Tr 



'Pi 
2m 



(A2) 



From the Lehmann representation we obtain the matrix element of gr~gaii^ some (arbitrary) A-representation: 

<A|gf - g!\^') = -2ni 5{E - E„){A\n){n\A'), En = {n\H\n), (A3) 



where |n) are the exact eigenstates of the unaveraged Hamiltonian. In a (spinless) disordered system, the states \n) 
are non-degenerate due to the disorder; thus one can rewrite llA3t in the |e)-basis with \n) = [e): 



Tr 



- S')] = \ E<^l£l^><^l^' - = E - e'){e\t-\e){e\e'W\e) = 6(E - e'){e'\^\e'h 



where d is the dimension (d = 2 in 2D). Substituting l|A4|l into (|AHl we obtain 



nco) = a{co)^Al = Tr[6«] = J dEv/^<E||-|E), a, ^ 0. 



Thus the main (Drude) contribution to the conductivity per spin projection is given by 



Od = - 



2e^ 
dmH 



J dEv/^<E||-|E). 



(A4) 



(A5) 



(A6) 



The easiest case is equilibrium at T = 0: then = -6(E) and the integral is equal to ^(Ol^j^jlO) s vp, so that 
od = pi ■ e^/h ■ inv/dm. 

Consider a wire of length L between two leads imder the voltage V, so that the energy distribution in the leads is 
given by f^''". If the effect of the interaction between electrons in the wire is weak, its distribution function linearly 
depends on the coordinate 



.20 



(AT) 



From l|A6t we conclude that also the leading (Drude) contribution to the conductivity is slightly inhomogeneous. 
Consequently, the stationary distribution of the electric field in the wire will also be inhomogeneous, so that the wire 
will be homogeneously charged. Thus the charge-neutralit y of the current-currying wire is slightly violated. 



APPENDIX B: GENERATING DIAGRAMS ON 
COMPUTER 

The diagrams are generated in the program^^ from the 
Kubo formula ll37l l according to the following algorithm. 



1. The simplest diagram is obtained from l|37l l by sub- 
stituting GF-operators (G^ and G^""') with their 
averaged values. The result corresponds to the 



diagram in Fig. 2(a) 



2. Add one cooperon or one diffuson line to the dia- 
gram(s) obtained on the previous step in all possi- 
ble ways. 

3. Leave only diagrams having no more than two 
loops (to calculate the number of loops, the dia- 



gram has to be redrawn in the "coordinate repre- 
sentation", see Sec. I VI Al l. 

4. Recursively perform two previous steps, until the 
last step produces no more new diagrams. 

5. Redraw all diagrams in the "coordinate represen- 
tation". Add single impurity lines in the Hikami 
boxes (in all possible ways), cf. dTOt . 



E.g., consider the diagram in Fig. |4(a)| without 
lonely disorder lines inserted in the Hikami boxes; 
each of its three diffuson lines is given by the in- 
finite series (l48l l. Leaving only the first terms in 
these series, we note that the resulting diagram is 
exactly the one depicted in Fig. |l(a)| with three dis- 
order (dashed) lines left. We see that the diagram in 
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Fig. |4(a)| includes the contribution which has been 
already taken into account in the weak localization 
diagram. To prevent double-counting, it should be 
subtracted. Such situation (when two different dia- 
grams have common contributions) can occur also 
for many other diagrams. All the diagrams must 
be checked for that; double-counted contributions 
must be subtracted. 

In principle, the above steps can be done manually, but 
it is better to make use of the computer program^*. On 
the first step we started with one diagram; in the end we 
get 215 onesi^ 



APPENDIX C: Sj^ AT co = 0: CALCULATING 
J Sk/{2nnfJ A^qKlnhf ON COMPUTER: ANALYTICS 
AND NUMERICS 



Once all diagrams with no more than two CD-loops 
are generated, they are automatically divided into differ- 
ent groups according to the number of HBs and diffuson 
lines (in total we get eight groups). The Hikami boxes 
of all diagrams are automatically calculated. We have 
checked that the diagrams of only one (the fifth) group 
contribute to a^n- The diagrams of the fifth group are 
calculated in the directory 2D . 5 / . 

Within the same group, the Hikami boxes can be 
summed up; after that the expression for the sum of all 
diagrams from the group can be written in the form l (80|l , 
with H'/f; being the sum of all HBs. 

Then we integrate over the diffuson momenta. We use 
dimensionless momentum variables (K and Q) defined 
in | [60l l. We checkedi^ that only the diagrams in Fig. |4] 
contribute to Sjg. Each of these diagrams has three diffu- 
son/cooperon lines; we label their momenta as K, Q, and 
K + Q. In case of small anisotropy of the energy spec- 
trum, 6 <^c 1, the denominator of the diffuson/cooperon 
depends only on the modulus of its momentum. Thus 
the angular dependence of the denominator comes only 
from (K + Q)^ - K^ + Q^ + 2KQ cos ip, where ip is the angle 
between K and Q. 

WecalculateJ d^k/{2nfi)^ J d^ t?/ (2 Tifi)^ in polar coordi- 
nates. In total there are two angular integrations; since 
the denominator depends only on one angle, another 
angular integration can be easily (and analytically) per- 
formed. The second angular integration is more com- 
plicated: the integrand is given by the sum of rational 
functions which have the form 



pin 

Jo 



di/) Pi(sini/', cos ip) 
2n P2(cos \p) 



(CI) 



where Pi 2 are polynomials. The denominator is even in 
ip, so we can leave only even (in i/') part of the numera- 
tor. Then numerator thus can be expressed as a another 



polynomial: 

/ 

Jo 



dipPaicosip) M V n , 

TT -51 T\' ^3 cos 4') = / . a» cos" \p. 

2n P2{cosu>) t—i 

^ ^' n>l 



Next, we perform the analytical integration over ip. Be- 
cause of the large number of terms to integrate and large 
size of the expressions, this analytical integration (basi- 
cally, calculation of residues) can only be done on com- 
puter. We use the fact that the denominators of all mass- 
full elements of the diffuson D^+q can be factorized into 
two expressions. In the programi^ ^^^7 denoted as 
uno= 1 + (K + Qf and due= 2 - (K + Q)^ + (K + Q)". The 
size of the integration result grows rapidly with powers 
of uno and due in the denominator, so it is necessary to 
split integrands into elementary fractions. ^-S' 

We assume that integrals over (k, q) converge on the 
scale oik < xH/l and q < xH/l, so that the integrals in 
dimensionless variables (K, Q) converge on the scale of 
K < and Q < 1. Together with the assumption x <^ 1, 
this permits us to use Taylor expansions (e.g., for Hikami 
boxes) in powers of W/fi 1 and ql/H 1. The usage 
of Taylor expansion here corresponds to the diffusive 
limit and is justified by the fact that we do not get any 
large-momenta divergences. 

Next, the integrand is symmetric with respect to 
K <r-> Q, and the integration operator has this symme- 
try too. So we symmetrize every term of the integrand, 
and express it in terms of new variables P = + 
and A - 2KQ/P. Accordingly, our integration operator 
is changed: 



dK dQKQ = 2 dK dQKQ 
Jo Jo Jo 



(C2) 



and / is the Jacobian: 



/ = 



d{P,A) 













= 4 Vl -A2. (C3) 



The integration results produced by integrate .max for 
A — > have insufficient precision (since the difference 
of large numbers has low numerical accuracy), so we 
integrate by P and A not from zero but from 0.001. Since 
VP the integrand^ at A = 0, and maximal contribution 
to the result occurs for A — > 1, this adjustment of the 
lower limit introduces only a negligible error into the 
result.— 

We treat separately (see files KpQ . max) the terms con- 
taining massless elements of Dk+q, and see that no diver- 
gences at small (P, A) occur, as it is predicted according 
to the VW- theorem in Appendix iDl 

The Jp dP together with the subsequent dA is per- 
formed in directory 2D . 5/A/. (As a check, we performed 
the numerical integration using alternative variables in 
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value at A^l: 0.334 




FIG. 7: The integrand for C 



the directory 2D . 5/B/, which led to the same result.) We 

/•oo pi 

use functions qagi for dP and qaws for dA, which 
are part of the quadpack package^ In order to minimize 
possible rounding corrections, all calculations are done 
with high precision (35 digits). 

The integrand for is plotted in Fig.[Zl The final 

result l(83]| has been calculated using different lisp real- 
izations (clisp,gcl and sb cl ), and for different numbers 
of digits in the numerical integration (20, 25, and 35). 

In conclusion, the steps of the diffuson momentum 
integration are: 

1. Perform the first (easy) angular integration. 

2. Select terms that contain massless component 
^k+q' Iritegrate them by (second angular inte- 
gration variable). 

3. Split other terms into elementary fractions, in order 
to decrease powers in their denominators. Inte- 
grate them by ip, and combine them with the terms 
obtained in the previous step. 

4. Calculate the integrand in many points within the 
interval A e [0,1]. Interpolate (linearly) the in- 
tegrand in the integration interval A e [0, 1], and 
integrate the interpolation result by A. 



Then, the Vollhardt-Wolfle theorem^ can be generalized 
to spin-dependent Hamiltonians: 

Theorem 1 In the the calculation of linear response coeffi- 
cients, no diffuson-type divergences occur if the unperturbed 
system possesses time reversal symmetry. 

In particular, diffuson singularities cannot occur in the 
calculation of the conductivity or spin susceptibility, as 
well as in the calculation of corresponding cumulants. 

We consider an arbitrarily complicated diagram with 
one or more diffusons. We choose any of them and prove 
that the coefficient in front of D^^ vanishes when the mo- 
mentum flowing through this diffuson approaches zero. 
Let us draw only the selected diffuson in the "coordi- 
nate representation" (see Sec. I VI At : other diffusons and 
all cooperons (if present) remain in the "ladder represen- 
tation"; e.g.. Fig. [S] Drawn in this way, a diagram con- 
sists of two bubbles with a wavy diffuson line between 
them [e.g.. Fig. 8(a)| . In a diagram for linear response 
to applied electric field, (at least) one of two vertices 
(we assume that it is the rhs-one) is proportional to the 
renormalized velocity operator lf59] |. Let us downgrade 
all diffusons/cooperons to crosses [removing disorder- 
averaging (dashed) lines connecting them] 42 After that 
the rhs-bubble of an arbitrary diagram can be written as 

Gr(p-q.p) 



(Dl) 



G!;"(p,p-q) 
d2p 



where Gi"''(p 



R vr 1'p) ^ri'i G^'(P'P ~ l) unaveraged 
Green functions in the with and «th order of the pertur- 
bation theory in the disorder potential. [Note that we 
did not draw the crosses on the GFs-lines in dPlb .l Since 
the disorder cannot break the time reversal symmetry of 
the imperturbed Hamiltonian, 



APPENDIX D: GENERALIZATION OF THE 
VOLLHARDT-WOLFLE THEOREM FOR 
HAMILTONIANS WITH SPIN-ORBIT INTERACTION 

We assume that the following statements (demon- 
strated above for the case of Rashba and Dresselhaus 
SOI) are valid in general for a system with time reversal 
symmetry: 

• The anomalous part of the velocity operator van- 
ishes at 6J = due to the vertex renormaliza- 
tion I 



The onli^ divergent (massless) element of the dif- 
fuson matrix ||49]| is D"". 



'im, n, pi, p2 



02 [G*^"Vp2, -pl)]^ 02 =G*f ^(pl, p2), 
Ol [G^'Vp2, -Pl)]^ 02 =G^'(pi, P2), 



(D2) 



which is a generalization of ll52t . From our assumptions 
it follows that Gr/a commutes with the renormalized 
velocity vertex vA. Applying the transformation l|D2|l to 
GFs in l|Dlt , transposing matrices under the trace, and 
substituting p -p in the integral we obtain 



-/ (S^P^sJS.^''^i^"(P-i'P)^A(p^p-q)] 
= /(^P^sJL^^^K'(P'P-'i)^A(p-q^p)l^ 



(D3) 
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FIG. 8: Illustration to the proof of the Theorem [T] Only one 
diffuson (with momentum q) is depicted in the "coordinate 
representation"; this diffuson Dq divides both diagrams in two 
parts (bubbles). The vertex on the rhs of both diagrams is 
proportional to the renormalized velocity operator {59\ . The 
rhs-bubble of (b) is the mirrored (and then rotated by 180°) 
rhs-bubble of (a). The sum of two diagrams is regular at £j — > 0. 

so that dPll l vanishes at q - 0, and the diffuson diver- 
gence is regularized. The proof can be performed for ev- 
ery diffuson in the diagram, so that all diffusons are reg- 
ularized. An example of mutually cancelling diffuson- 



div ergen ces is depicted in Fig. |8l both diagrams |8(a)| 
and |8(b)| diverge at q — > 0, but their sum is regular. A 
similar cancellation takes place for the diagrams in Fig.|3l 
the singularity in D'^ in the diagr am |4(a)| is cancelled (i.e., 
regularized) by the diag ram |4(b)[ while the singular ity in 
Dq" in the diagram |4(a)| is cancelled by the diagram |4(c)[ 
Such cancellation is provided by the fact that in a sys- 
tem with time-reversal invariance cooperon components 
are exactly the same as components of a diffuson. (This 
remains true also in the presence of interaction between 
electrons^!) 

Thus we demonstrated that (in the absence of inter- 
action) there can be no diffuson type singularities in the 
presence of the time-reversal symmetry, thus generaliz- 
ing the theorem proved in Ref . |^ for the spinless case 
without SOI. 
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Symbol 


Its definition 


d 


spatial dimension 


yi 




CandD 


see JID and JIS 


Gr and Ga 


see Gl; = G], 


I 


T n r 1 — r 1 — tz. 

mean free path oi an electron between 

\D lAL/OvT^J U.tr± L L tritU L±^ J\w.Cl L Ltr± li. Ill ±±± I vj LAX X LXt^ J 




electron (orbital) dephasing length due to 
inelastic scattering 




dephasing rate due to inelastic scattering 


f 


(temperature-dependent) chemical potential 


Pf = -\j2mj.i 


Fermi momentum 


Vf = pf/m 


Fermi velocity 


V, DoS 


density of states 


1Rand3 


real and imaginary part 


ais and aan 


see and 00) 


S 


see IZD 




2x2 unity matrix 


01, 02, and 03 


Pauli matrices 




see {ID 


Pf , Xa, and X,, 


see fist 


X and b 


see 113 


CD 


cooperon or diffuson 


CS 


coordinate system 


DM 


density matrix 


GF 


Green function 


GFB 


Green functions box 


HB 


Hikami box 


Ihs 


left hand side 


rhs 


right hand side 


SOI 


spin-orbit interaction 


vw 


VoUhardt-Wolfle 


WL 


weak localization 


ZLA 


zero-loop approximation 



TABLE I: Some often used notations and abbreviations. 



